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Abstract. - We use confocal microscopy to study the motion of a magnetic bead in a dense 
colloidal suspension, near the colloidal glass transition volume fraction (j>g. For dense liquid-like 
samples near (f)g, below a threshold force the magnetic bead exhibits only localized caged mo- 
tion. Above this force, the bead is pulled with a fluctuating velocity. The relationship between 
force and velocity becomes increasingly nonlinear as (f>g is approached. The threshold force and 
nonlinear drag force vary strongly with the volume fraction, while the velocity fluctuations do 
not change near the transition. 



Introduction. - Near the glass transition, molecular motion is greatly slowed, yet the 
motion of molecules is difficult to directly observe and the character of the molecular slowing 
is hard to determine [1-4]. Thus, colloidal suspensions of micrometer-sized spheres are a 
useful model system for investigating the nature of the glass transition, as they can be directly 
observed using optical microscopy [5-9]. As the volume fraction of a colloidal suspension is 
increased, the motion of the particles becomes restricted due to the confining effect of the 
other particles: each particle is temporarily trapped in a "cage" formed by its neighbors 
[9-12]. For liquid-like samples, these cages eventually deform and allow particles to diffuse 
through the sample [9, 10]. Above a critical volume fraction {(f)g ~ 0.58 for hard sphere 
colloids), the particles' diffusion constant becomes zero and the macroscopic viscosity increases 
dramatically [10,13]. A system at such state is considered a colloidal glass and 0g is identified 
as the colloidal glass transition point. To gain insight into the nature of the colloidal glass 
transition, past work used optical microscopy to study the motion of individual particles as 
a function of volume fraction [5-9]. These experiments discovered that particle motion in 
equilibrated samples is both spatially heterogeneous and temporally intermittent. However, 
only a few studies [14, 15] have focused on nonequilibrium, driven motion. 

In this Letter, we study experimentally the motion of a magnetic bead pulled through a 
dense colloidal suspension. Similar to previous work with magnetic probes, we are able to 
probe the local rheological properties of the suspension [16-18], although unlike those studies, 
our samples are heterogeneous on the scale of the magnetic probe. In our experiments, as the 
colloidal glass transition is approached the motion of the magnetic bead becomes complex. 
Below a threshold force, the magnetic bead exhibits only localized caged motion. When the 
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external force is above this threshold, the magnetic bead moves significantly faster, with a 
fluctuating velocity. Closer to 4>g, the threshold force rises, and the velocity- force relationship 
becomes increasingly nonlinear. Our measurements reveal dramatic changes in the drag force 
acting on the magnetic bead as the glass transition is approached. 

Experimental. - The colloidal suspensions are made of poly-(methylmethacrylate) parti- 
cles, sterically stabilized by a thin layer of poly-12-hydroxystearic acid [19]. The particles have 
a radius a = 1.10 ^m, a polydispersity of '^5%, and are dyed with rhodamine. The particles 
are slightly charged, shifting the phase transition boundaries from those of hard spheres. We 
observe the freezing transition at 0.38, the melting transition at (pm ~ 0.42, and the glass 
transition at (j)g ~ 0.58. The glass transition is characterized by a vanishing diffusion constant 
for the particles' motion, measured by confocal microscopy. The colloidal particles are sus- 
pended in a mixture of cyclohexylbromide/ cis- and trans- decalin which nearly matches both 
the density and the index of refraction of the colloidal particles. We add superparamagnetic 
beads with a radius of omb — 2.25 /im (M450, coated with glycidyl ether reactive groups, 
Dynal), at a volume fraction 10^^, so that the magnetic beads are well separated. We do 
not observe attraction or repulsion between the colloidal particles and the magnetic beads, in 
either dilute or concentrated samples. Before beginning experiments, we mix the sample with 
an air bubble which breaks up any pre-existing crystalline regions, then we wait 20 minutes 
before taking data, to allow the air bubble to stop moving. The colloidal suspension is not 
observed to flow during the acquisition of the data. 

With a confocal microscope [20], we rapidly acquire images (1.88 image/s) of area 80 
/um X 80 /im, containing several hundred particles. The magnetic beads are not fluorescent 
and thus appear black on the background of dyed colloidal particles. To exert a force on 
these beads, a Neodymium magnet is mounted on a micrometer held just above the sample. 
The force is calibrated by determining the velocity of the magnetic beads in glycerol for a 
given magnet position, and inferring the drag force from Stokes' Law {F = QirrjaMBV, where 
F is the drag force, rj is the viscosity, and v is the observed velocity) [17]. The imperfect 
reproducibility of the magnet position causes a 5% uncertainty of the force. Additionally, 
the variability in magnetic bead composition results in an uncertainty in the magnetic force 
on different beads, which we measure to be less than 10%. Also, the magnetic beads are 
not density matched, and their effective weight is 0.1 pN. This is nearly negligible compared 
to the applied horizontal forces in our experiments, and so in our data below we consider 
only the applied forces and the measured velocities in the x direction. We observe very little 
vertical motion of the magnetic beads, typically less than 5% of the horizontal motion. This 
is less than our other measurement errors, discussed below. We study isolated magnetic beads 
at least 35 /im and more typically >50 /im from the sample chamber boundary and from 
other magnetic beads. Our experiments are limited by crystallization as the average velocity 
of the magnetic bead decreases significantly in the crystallized regions. Thus, we make sure 
that the data are collected well before crystallization appears. Repeated measurements are 
reproducible before the onset of crystallization. 

Results and discussion. - We study the motion of the magnetic beads for a range of 
forces. Within each image, we locate the center of the magnetic bead with an accuracy of at 
least 0.5 /im [21]. To show the behavior of the magnetic bead, we plot the average magnetic 
bead velocity against the applied force on a log-log plot in Fig. ^ for a range of volume 
fractions. The average velocity of the magnetic bead for a given F is given hy v = Ax/ At for 
the entire course of the measurement; the error bars in Fig. ^ are from the standard deviation 
of repeated measurements of v for the same magnetic bead at the same F. For a dilute system 
{(j) = 0.29, asterisks) the data fall on a hue with a slope of 1, indicating that the velocity is 
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Fig. 2 



Fig. 1 - Average velocity of the magnetic bead in the direction of the motion as a function of force 
acting on the magnetic bead for the volume fractions: 0.29 (*), 0.45 (A), 0.50 (t>), 0.52 (o), 0.53 (o), 
0.55 (<). Solid lines are fits to the data obtained using Eq. (1). Each data set is taken for a single 
magnetic particle and thus the magnetic particle variability may shift curves relative to each other 
but does not influence the shape of a curve for a given volume fraction. The velocity error due to the 
uncertainty in the bead position is 5%, smaller than the symbol size. The slope=3 line is added for 
comparison. Inset: Data for </> = 0.29 on a linear/linear plot. 

Fig. 2 - Fitting parameters as a function of volume fraction: (a) Vb, (b) threshold force Fq, and (c) 
power exponent a. The error bars are due to the uncertainties in measuring the forces and velocities 
that the model is fit to. The solid symbols for Fq correspond to data sets for which stalled behavior 
was directly observed below Fo (all data sets (j> > 0.50). 



linearly proportional to the applied force and Stokes' law applies (Fig. ^ inset). Due to the 
nonzero volume fraction, the effective viscosity is modified from that of the solvent (770 = 2.25 
mPa-s) to the larger value of 77 = 6 ± 1 mPa-s, in reasonable agreement vifith past macroscopic 
measurements [13]. 

At higher volume fractions the velocity-force relationship becomes more complex (Fig. ^ . 
For data at large forces, the average velocity seems to grow as a power law with force, v ~ F", 
with exponent a between 1.5 and 3. From left to right the curves increase in volume fraction, 
and the slope at the highest forces for each curve increases nearly monotonically, reflecting 
the approach of the glass transition. For all samples (j> > 0.50, we took data at extremely low 
forces, and found that the magnetic bead did not move more than our resolution limit within 
four hours; instead, it exhibits random caged motion. This observation time is limited by the 
onset of crystallization of the sample; in other words, the velocity is smaller than 4- 10~^ /xm/s. 
Since the colloidal particles move due to Brownian motion one would expect rearrangements 
to "allow" for the magnetic bead to move, even if this motion would be primarily random with 
a small bias from the external force. However, within our observation time scale we cannot 
see this behavior, and conclude that the velocity of such a random walk is quite small. For 
(f) = 0.29 and 0.45, for all forces checked, motion of the magnetic bead was observed. (Forces 
below 0.1 pN are unobtainable, as this is the effective weight of the magnetic bead.) 

To determine how the velocity /force relationship changes as the volume fraction is in- 
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creased, we fit the data to the following equation: 

v= vo-{^-ir, F>Fo (1) 
0, F<Fo 

where Fq is a threshold force and a is the exponent characterizing the growth of the velocity 
at large forces, vq is the velocity the magnetic bead would have if the applied force F was 
equal to 2Fo . This functional form was also used to fit data from a two-dimensional simulation 
similar to our experiments [22]. The solid lines in Fig. ^ are the results of fitting the data 
with Eq. (1), with the exception of the 4> = 0.29 data (asterisks) where Stokes' Law applies. 
While the model seems a good fit to our data, our data are not strong enough to verify the 
model. However, the model provides a useful way to quantify the main features of the data. 

As the colloidal glass transition is approached the system becomes more viscous, and this 
is reflected in the changing fit parameters from Eq. (1) shown in Fig. |2] For low volume 
fractions, the three parameters wq, Fq, and a are all small, and Stokes' Law may be the 
simplest explanation for the force/ velocity data. The fit parameters change as the volume 
fraction (p (j)g. vq initially grows with increasing (j) and then drops [Fig.|3[a)]. The initial 
increase presumably reflects the change from the vq —^ limit at low volume fractions. The 
decrease as the glass transition is approached may be due to the sharply increasing effective 
viscosity acting on the magnetic bead, vq is also sensitive to variations in the other fitting 
parameters, as for large F, v ^ (woFg~")F". 

The threshold force Fq grows as increases, reaching a plateau value for large 0. We 
expect that at (j)g the threshold force will be finite, both from our data, and also as it seems 
more likely that Fq diverges at (/)rcp ~ 0.64, the random close-packing volume fraction [13]. 
The value of the threshold force near Fq ~ 1 pN, is similar to that predicted by the work of 
Schweizer and Saltzman et al. [11,12]. They determined an effective free energy characterizing 
caged particle motion, and the derivative of this energy can be used to find an effective force 
to break out of a cage, resulting in a value of w SOfc^T/a w 0.2 pN for (j) = 0.6 [11, 12]. This 
prediction considers the caged particle to be the same size as the surrounding particles, and 
for our larger magnetic beads the stall force would likely be larger, closer to our experimental 
values. Thus, our measured threshold forces Fq may be related to the strength of colloidal 
cages; previously only the sizes of these cages have been measured [9, 10]. 

The exponent a changes dramatically as the glass transition is approached, as seen in 
Fig. 121c), where it rises to almost 3 near (pg. Our results are in contrast to Ref. [22]; while the 
model fits their data as well, they found a — 1.5 for glassy systems. The difference may be 
due to the much different particle interactions: they studied unscreened point charges in a 2D 
system, qualitatively different from our colloidal particles. In our experiments, it is unclear 
if a diverges at (j)g or continues to grow as 0rcp is approached. Experiments at cj) > cjjg are 
difficult to interpret due to the aging of the sample [23] : for a given force the magnetic bead 
response depends on the age of the sample. 

The motion of the magnetic bead is not smooth, as is seen by a typical plot of its position 
as a function of time [Fig.|3Ia,b)]. While the largest "jumps" seen in Fig.|2fa) are roughly the 
diameter of the surrounding colloidal particles (2a ~ 2 fim) , at other times the magnetic bead 
takes much smaller steps. When the magnetic bead moves faster the fluctuations become less 
significant and the motion of the magnetic bead is smoother, shown in Fig. |2Ib). 

We arbitrarily choose At = umb/v as a time scale and define the instantaneous velocity 
as v{t) — [r(t + At) — r(i)]/Ai. This time scale is indicated with the scale bars shown in 
Fig. inia,b), and the results that follow are not sensitive to this choice (and in particular do 
not change using a rather than omb)- The instantaneous velocity is shown in Fig. |2Ic). The 
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Fig. 3 - (a) The position of the magnetic bead in the direction of the motion as a function of time, 
F = 0.58 pN, (/> = 0.52, v = 0.0041 fim/s. The scale bar corresponds to umb/v = 550s. (b) The 
position of the magnetic bead in the direction of the motion as a function of time, F = 6.46 pN, 
(j} = 0.55, V = 0.80 /im/s. The scale bar corresponds to umb/v = 3 s. (c) The instantaneous velocity 
in the direction of motion Vx and perpendicular to the motion Vy as a function of time. The velocities 
are calculated from the displacement data using At = 3 s. The dashed lines denote the average 
velocities Vx ~ 0.80 fim/s and Vy = 0.01 /im/s. The data are the same as (b). (d) The probability 
distribution of the instantaneous velocities Vx (o) and Vy (A). The solid lines denote a Gaussian fit 
with Gvx = 0.24 fim/s and a^y — 0.10 /im/s. The data are the same as (b). 



magnetic bead velocity has much larger fluctuations along the direction of motion [vx, top 
trace in Fig.|3Jc)], while fluctuations in the transverse direction are smaller. The distributions 
of Vx and Vy, shown on Fig. |3{d), are Gaussian (solid fit lines). The fluctuations in x and 
y appear to be only weakly correlated. The velocity fluctuations reflect the fact that the 
magnetic bead is of similar size to the colloidal particles: its motion is sensitive to colloidal 
particle configurations. Furthermore, dense colloidal suspensions are spatially heterogeneous 
systems [5,6] and some regions may be "glassier" and thus harder to move through. 

As the glass transition is approached, the spatial heterogeneity increases [5,6], and one 
might expect that the velocity fluctuations would become more significant. However, this 
is not the case, as is seen in Fig. 01 which shows the standard deviation of the instantaneous 
velocity ayx plotted against the average velocity v. Different symbols indicate different volume 
fractions, showing that it is not the volume fraction which determines the standard deviation 
but the average velocity: the fluctuations are larger when the average velocity is larger. Denser 
suspensions have slower velocities for the same force (as seen in Fig.^ but also correspondingly 
smaller velocity fluctuations. In Fig. ^ the data are more scattered at low velocities and for 
any given sample the data scatter on both sides of the fit line. For different choices of At, we 
find that ayx ^ with the exponent (3 = 0.8 ±0.1. The smaller the At the wider distribution 
of a^x- However, the general dependence of a^x on v does not depend on At. Wc find that 
avy ~ QAgvx and thus a similar trend holds for the transverse fiuctuations. 

The changing character of the motion as the glass transition is approached can be un- 
derstood by calculating the modified Peclet number. This is the ratio of two time scales, tb 
characterizing the unforced motion of the colloidal particles, and tm characterizing the mo- 
tion of the magnetic bead. The Brownian time tb — /2Dao is the time it takes for colloidal 
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V (um/s) 

Fig. 4 - The standard deviation of the instantaneous distribution of Vx, {(Jvx) as a function of «. The 
solid line is a fit to the data and has a slope of 0.80. The symbols represent the different volume 
fractions which are the same as Fig. Q 

particles to diffuse a distance equal to their radius a and is based on their asymptotic diffusion 
constant Doo- Doo becomes small near the glass transition and becomes zero for a colloidal 
glass. The magnetic bead time scale tm = a/v is the time scale for the magnetic bead to 
move the distance a and is based on the magnetic bead's average velocity v. The modified 
Peclet number is Pe* = tb/tm and is larger than 1 for all symbols shown in Fig.Q] [24]. For 
the largest forces, Pe* « f 0000. Thus, the forced motion of the magnetic bead is much more 
significant than the Brownian motion of the surrounding colloidal particles: the magnetic 
bead pushes these particles out of the way, plastically rearranging the sample. This results 
in a lowering of the effective viscosity acting on the magnetic bead at higher velocities. On 
the macroscopic scale, this may correspond to traditional rheological measurements at high 
strains. Past work found for similar colloidal samples that the viscoelastic moduli decreased 
as the maximum strain increased [25]. 

For extremely low forces, Pe < 1 and the velocity should be linearly related to the force. 
However, the velocities predicted by linear response theory are well below what we observe 
(and for low forces, below what we can measure). This is not surprising, as the lowest forces we 
can apply are fairly large; a nondimensional way to characterize this is the ratio ausF /ksT 
which is ~ 50 for F=0.1 pN. 

Since Pe* > 1 in our experiments, the Brownian motion should be unimportant, similar to 
a granular system. Moreover, conjectures of the existence of a "jamming transition" speculate 
that the colloidal glass transition is similar to jamming in granular media [26], and so it 
is interesting to compare our results with studies performed in a granular media. Albert 
et al. [27] immersed a cylinder in granular particles and measured the force exerted on the 
cylinder when it moved with constant speed relative to the particles. They found that the 
drag force on the cylinder is independent of the velocity, quite different from our behavior 
(Fig. In our experiments, the solvent viscosity may be important, causing the velocity 
dependent drag force. 

In conclusion, we find that near the glass transition the motion of a microsphere through 
a colloidal suspension changes dramatically. A threshold force appears, below which the 
motion of the sphere is localized, and this threshold force increases as the glass transition is 
approached. The existence of the threshold force hints that the system may locally "jam" 
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even when the colloidal suspension is globally still a liquid. Above the threshold force the 
velocity is related to the force by a power law and the sphere moves with a fluctuating velocity, 
locally deforming the sample. At high forces, the velocity is related to the force by a power 
law. This can be inverted to obtain an effective drag force on the sphere, growing weakly with 
velocity, F ~ v^/" with a growing from 1 far from the glass transition to nearly 3 close to the 
transition. These results indicate that the approach of the colloidal glass transition is signaled 
by a growing nonlinearity of the drag force. This may be a new way to characterize the glass 
transition and should serve as a useful test for some glass transition theories [11. 12.28,29]. 

We thank R. E. Courtland, M. B. Hastings, S. A. Koehler, D. Nelson, K. S. Schweizer, 
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